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Abstract 

Adhesion is a key issue for researchers of various 
fields, it is therefore of uppermost importance to 
understand the parameters that are involved. Com- 
monly, only surface parameters are employed to 
determine the adhesive forces between materials. 
Yet, van der Waals forces act not only between 
atoms in the vicinity of the surface, but also be- 
tween atoms in the bulk material. In this review, 
we describe the principles of van der Waals inter- 
actions and outline experimental and theoretical 
studies investigating the influence of the subsur- 
face material on adhesion. In addition, we present 
a collection of data indicating that silicon wafers 
with native oxide layers are a good model substrate 
to study van der Waals interactions with coated 
materials. 

Introduction 

Stimulating or preventing adhesion is a key issue 
for researchers of various fields. To solve these 
problems, a comprehensive understanding of the 
prevailing adhesion mechanism is indispensable. 
Yet, not only various adhesions mechanisms ex- 
ist, but also plenty of parameters that can affect 
adhesion: Nanoscale or microscale roughness,'^ 
static charges or the zeta-potential at the inter- 
face,^^ surface energies,^ ^ and contact shapes'^ 
are a few frequently-studied examples. All these 
parameters, however, have in common that they 
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are describing the surface of a material. Hence, 
the question arises whether adhesion is really only 
'superficial'. This question is of great importance 
since commonly used photoresists, coatings, ad- 
hesion promoters or other functionalized surface 
layers are often in the range of just a few nanome- 
ters. These dimensions are smaller than the range 
of interactions such as van der Waals (vdW) in- 
teractions. Hence, the material right underneath 
the surface might indeed have an effect on adhe- 
sion mediated by vdW interactions. In this paper, 
we review experimental and theoretical studies in- 
vestigating the influence of subsurface material on 
vdW forces. Additionally, we provide a collec- 
tion of experimental data highlighting the suitabil- 
ity of stratified substrates based on silicon wafers 
to study vdW interactions. 

Van der Waals interactions 

Already more than hundred years ago van der 
Waals introduced his theory of an attraction be- 
tween neutral atoms in order to explain non-ideal 
gases.l^ Later, three types of interactions were 
identified to contribute to the vdW interactions: 

Keesom interactions characterize dipol-dipol in- 
teractions of molecules that carry permanent 
dipoles.'^ 

Debye interactions describe forces between a 
permanent dipole that induces a dipole mo- 
ment in an otherwise unpolar molecule.'^ 
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Figure 1: VdW potentials ^{d) for different geometries. Adapted from.'^ 



London interactions, also called dispersion in- 
teractions, describe forces between instanta- 
neously induced dipoles.'^ 

All three parts have in common that the interac- 
tion energy scales with — ^ . Hence, the vdW po- 
tentials ^ for the interactions between two single 
atoms separated by a distance d can be written as 



^{d) = -C/d^ 



(1) 



Due to this scaling, vdW interactions are often 
considered as of short-range. 

Hamaker, however, calculated energy-distance 
relations for macroscopic objects by pairwise sum- 
mation over all atoms, continuing the work of 
Bradley and DeBoer.!^^^^ Depending on the ge- 
ometry, different scaling laws apply (cf. Figure[T]). 
To account for the properties of the involved ma- 
terials, Hamaker introduced a coefficient A, also 



called 'Hamaker constant', which he defined to 

A = K^CpiP2 (2) 

where p,- are the number of atoms per unit vol- 
ume of the two materials. The controversial sub- 
ject of Hamaker's theory was that he assumes a 
pairwise additivity of the vdW interactions, which 
is generally speaking not valid. A few years 
later, Casimir used a completely different ansatz 
to calculate the force between two ideally conduct- 
ing semi-infinite half-spaces in vacuum.^ On the 
basis of Planck's famous theory, he summed up 
the allowed electromagnetic modes between two 
conducting plates. Lifshitz extended Casimir's 
idea and presented a theory for arbitrary materi- 
als, based on quantum field theory.!^ In princi- 
ple, although many studies differentiate between 
Lifshitz-vdW and Casimir interactions, Lifshitz 
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and Casimir essentiall y desc ribed the same effect, 
but with different foci.'^^^ By treating the inter- 
acting objects as continuous media, Casimir's and 
Lifshitz' theories circumvent the question of pair- 
wise additivity. Interestingly, Lifshitz' ansatz led 
to the same scaling laws as the classical Hamaker 
ansatz (cf. Figure [T]). Only the way the Hamaker 
constants are derived is different. Following Lif- 
shitz' theory, they are calculated from the optical 
properties of the involved materials and can be ap- 
proximated (see Appendix) by 
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with the dielectric constants £;, the refractive in- 
dices in the visible regime n,, and the main elec- 
tronic absorption frequency We. 

The works of Hamaker, Casimir and Lifshitz 
demonstrate that vdW-Casimir interactions can in- 
deed be regarded as of long-range, since for meso- 
scopic and macroscopic objects, the absolute value 
of the exponent of the scaling law is decreased (the 
interaction between two semi-infinite half slabs 
scales with — for instance). Yet, the long-range 
character is restricted due to the finite speed of 
light.'^ This retardation effect increases the abso- 
lute value of the exponent of the scaling law by up 
to one (for d ^ 20 nm). For separations smaller 
than lOnm, howe ver, the retardation can usually 
be neglected.E^ 

Van der Waals interactions with 
coated substrates 

Using the equations given by the theories men- 
tioned above, it is usually possible to predict the 
potentials for the interactions of two uniform ob- 
jects. Yet, many systems consist of coated sub- 
strates. However, VdW interactions act not only 
between atoms in the vicinity of the surface, but 
also between atoms in the bulk material. Early 
experiments of Israelachvili und Tabor, showed 



that in the case of mica substrates covered with 
a monomolecular layer of stearic acids, both the 
acids and the mica contribute to the vdW in- 
teractions.^^ These experimental studies were in 
agreement with theoretical predictions of Lang- 
bein, who postulated that the interactions with the 
surface layer dominate for separations d smaller 
than the layer thickness D (D <^ d) and the in- 
teractions with the bulk mate rial dominate in the 
opposing limit (D ^ |22 | 23 | yiQj-Q j-gcent studies 
showed that variations in the thickness of a sur- 
face layer induce differences in the vdW poten- 
tials and influence e.g. the stability of thin liquid 
coatings. ^^ ^^ These thin film dewetting studies 
moreover demonstrated quantitatively the impact 
of the subsurface composition to the effective in- 
terface potential:'^ The impact was measured ex- 
perimentally by determining the differences in the 
preferred wavelength of spinodally dewetti ng thi n 
films with variable subsurface composition.l25E2l 

Interactions in biological systems 

VdW interactions also play a major role in bio- 
logical systems. Particularly non specific adhe- 
sion is governed by vdW interactions in conjunc- 
tion with electric double layer inte racti ons, usually 
described using the DLVO-theory^'^^' or extended 
DLVO-theory.^^ Many biological processes, such 
as the aggregation of proteins, the unspecific ad- 
hesion of cells and bacteria,'^^^ the adherence of 
abalones,!^ and the sticking of geckosl^^^®* are in- 
fluenced and sometimes dominated by vdW inter- 
actions. As a logical consequence, these processes 
are also not 'purely superficial', viz. not solely de- 
pendent on the properties of the surface. It could 
be shown that, when interacting with a coated 
substrate, proteins indeed sense both, the surface 
layer and the underlying material: Adsorption ex- 
periments on tailored silicon wafers with differ- 
ences in the thickness of the oxide layer on top 
of the wafers revealed qualitatively different ad- 
sorption kinetics of multiple types of proteins. ^^^^ 
Based on Monte-Carlo simulations featuring sur- 
face processes such as surface mobility and con- 
formational changes, the distinctions were invoked 
by the influence of the vdW-interactions on the 
time scale of these processes. X-ray reflectivity 
experiments corroborated these findings. "^"^ Recent 
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Figure 2: Schemes of two multilayer configurations: A) symmetrical triple layers and B) interaction 
between a uniform material and a coated material. 



studies discovered that the influence of the sub- 
surface material on adhesion is sensed by larger 
biological objects, too: The unspecific adhesion 
of bacteria from the Staphylococcus genus is af- 
fected by the properties of the subsurface ma- 
terial, as could be shown by AFM force spec- 
troscopy measurements.^^ Moreover, adhesion ex- 
periments with setal arrays of live geckos revealed 
that the adhesion force was significantly varied 
by a change in the substrates subsurface compo- 
sition.!^ 

Theoretical description of multilayer 

systems 

For a comprehensive theoretical description of the 
vdW potentials for systems involving multilayer 
structures, not only the surface layer, but also the 
composition of the entire substrate has to be taken 
into account. Unfortunately, there is no general 
equation for the vdW interactions between arbi- 
trary objects and layered substrates. In the follow- 
ing we therefore will focus on the non-retarded in- 
teractions between infinite planar interfaces. 

Generally, a common way to find out unknown 
Hamaker constants is the use of combining rules 
(geometric mean).'^ These relations, that are de- 
rived from the combining rules for surface ener- 
gies, may be used to calculate 'effective' Hamaker 
constants for a multilayer system. Yet, these rela- 
tions break down whenever the Keesom and De- 
bye part (the zero frequency terms) cannot be ne- 
glected. ^ Especially in multilayer systems, where 
multiple Hamaker constants are necessary, com- 



bining rules are not applicable. On the basis of the 
Lifshitz theory, the potentials for the interactions 
of symmetrical triple layer films (cf . Figure [2]A) 
were calculated by 
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with Aij^ki the constants for the interactions of 
the two different interfaces.'^^^^ Using the same 
ansatz, the vdW potential of the interactions be- 
tween a probe material and a substrate coated with 
a layer of thickness D (cf. Figure [2]B) is given 
byM 
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For larger separations D <^ d, however, Eq. ([5]) 
is no longer valid. For the description of an exper- 
imental system with a variable D (e.g. a thickness 
of a coating), we hav e prev iously chosen an alter- 
native approximation :'25I261 gy assuming a scaling 
of the interaction with % -|- and a continu- 
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ous transition between the boundary cases D <^ d 
and D:^ d,we gained 
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where A 12-42 stands for the interaction of material 
1 via medium 2 with material 4 in the case that 
D — 0, viz. medium 3 is nonexistent. 
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Figure 3: Properties of thin oxide layers: A) Theoretically expected transition of the refractive index (short 
dashes) and a proposed double layer oxide model (long dashes). B) Single layer oxide models whereby i) 
the oxide layer has i) Si02 bulk-like properties (solid line) or ii) properties expected by optical studies with 
limited resolution (dotted line). C) Calculated potentials of the van der Waals interactions between a semi- 
infinite half-space of polystyrene (^vis = 1.585, £o = 2.6) and a silicon wafer with a thin oxide layer for 
the three model configurations (cf. A and B) and a bulk silicon dioxide substrate (Wyis = 1-46, £o = 3.^^ 
as a comparison (dash-dotted line). To include the short-range repulsive interaction a second term ^ 
(maintaining the difference in exponents of the Lennard- Jones potential) was added. The parameter C was 
kept constant for all configurations, since the surface is essentially the same. 



Silicon wafers as a model system 

Most of the studies mentioned above used sili- 
con wafers as a model system to study vdW inter- 
actions with coated materials. As described ear- 
lier, these interactions are essentially dependent 
on the optical properties of the involved materi- 
als (Eq. ([3])). The properties of silicon oxide films, 
however, have been discussed controversially for 
decades. In the following, the term "silicon oxide" 
stands for SiO^. The use of Si02 bulk values for 
oxide films thicker than 100 nm is generally ac- 
cepted. Yet, the validity of these values for thin 
film has been questioned numerous times and es- 
pecially the optical properties, such as the refrac- 
tive index, are disputed. In general, thin silicon ox- 
ide layers may be described by two different mod- 
els: 

Single layer model: By assuming a sharp transi- 
tion between the bulk Si and the oxide, the 
latter can be described by a single layer (cf. 



Figure|3]A). In this case, the oxide layer may 
have the same refractive index as bulk Si02 
or an increased refractive index. 

Double layer model: Since a sharp transition be- 
tween the materials is highly unphysical, a 
continuous transition or an interface rough- 
ness is very likely. As the thickness of 
this transition region is of the same order 
of magnitude as the oxide layer thickness, 
a double layer configuration is a obvious ap- 
proximation (cf. Figure [3]B). 

The van der Waals interactions sensed by an probe 
object should differ strongly for these different 
model configurations (cf.Figure[3]C). 

By applying the single layer model, Jellison 
observed an increase in the refractive index of the 
whole oxide layer. Using spectroscopic polar- 
ization modulation ellipsometry, the refractive in- 
dex of very thin layers was determined to 1.5 — 1.8 
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Table 1 : Surface properties of a native oxide layer and of a thermally grown thick oxide layer as a com- 
parison: thickness (d), root mean square (rms) roughness (by (l/im)^ AFM scan), advancing (adv) and 
receding (rec) water contact angle, and surface energy 7 (obtained from contact angle measurements of 
three different liquids ^2). 

oxide layer d [nm] rms [nm] 0adv n 0rec[°] 7 [mJ/m^] 

native 1-7(3) 0.13(3) 5(2) compi. wetting 63(1) 
thick 150(1) 0.09(2) 7(2) compi. wetting 64(1) 



(at A = 800 nm). These findi ngs were later on con- 
firmed by other studies. '^^'^ 

Experimental support for the double layer 
model, viz. the observation of an interfacial transi- 
tion layer, is also given by previous studies. Exper- 
imental studies noticed a thin (^ 6 — 7 A) region 
of atomically mixed Si and O with an refractive 
index of n ^ 2.8 - 3.2 (at 546.1 nm).i22l58| High- 
resolution core-level and XPS spectroscopylSSH 
also confirmed that "the interface is not abrupt, 
as evidenced by the high density of intermediate- 
oxidation states (about two monolayers of Si) and 
by their nonideal distribution" .^SH These findings 

were matched by predictions of theoretical mod- 
els. ^2 63 

The results of theses studies, however, are not 
contradictory, but arise from the different methods 
applied. A problem of optical reflectivity methods, 
such as ellipsometry, is that they are not able to de- 
termine the density and the thickness of thin films 
(< 5nm) independently (not to mention to distin- 
guish between two of such films). Thus, the usage 
of the single layer model for these methods is the 
only possible way. Yet, for thinner silicon oxide 
films, the transition layer fraction of the total ox- 
ide layer is increased resulting in an observation of 
a higher overall refractive index. 

Another limitation of all of the mentioned stud- 
ies is, that they are performed on sihcon wafers 
with thermally grown silicon oxide layers. Espe- 
cially for thin (< 10 nm) films, the process param- 
eters of the artificial growing process can have a 
significant influence on the density and the optical 
properties of the silicon oxide (decreased vs. in- 
creased refractive index). Yet, already without 
any pretreatment, silicon wafers are covered with 
a native oxide layer of 1.5 nm to 2nm thickness. 
Since only limited data is available for the optical 
properties of native oxide layers, we present a brief 
summary of the properties of these layers contain- 
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Figure 4: Electron density profiles of a native 
(dashed line) and a thick oxide layer (solid line) 
determined by high energy X-ray reflectivity. The 
z-coordinate has been set to zero at the center of 
the transition between Si and its oxide. 



ing previously unpublished data. 

Properties of thin native oxide layers 

The increase in refractive index to values up to n = 
1.8, as predicted by several studies, ^^^^ is highly 
unlikely for native oxide layers as polymer dewet- 
ting studies have shown via an indirect wa.y^ 
Thin liquid polystyrene (PS) films prepared on Si 
wafers with native oxide layers (Z)sio = 2.4nm) 
were unstable and dewetted spinodally. Since this 
process is driven by the minimization of th e free 
energy determined by the vdW potential'^'^the re- 
fractive indexes of the oxide cannot be higher than 
the one of PS 1.59), as shown by Eq. ([3]). 

A more direct ansatz is to compare native ox- 
ide layers to thick oxide layers in terms of the ma- 
terial properties, such as surface roughness, sur- 
face chemistry, homogeneity, electron density and 
stoichiometry: The surface characterization via 
atomic force microscopy (AFM) and contact angle 
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binding energy [eV] 

Figure 5: A) XPS core level spectra for Si-2p (Al-Ka radiation, hco = 1486.6 eV, normal emission) for 
different stages of surface ablation by Ar ion etching with asymmetric decrease of intensity contributions 
from oxide species. B) Double layer oxide model with Si*^ - (Si^+, Si^^, Si^+) - Si^^ stacking (from 
bulk to surface). C) Simulation of spectra for the double layer model from B) with the same asymmetric 
decrease of oxide intensities as observed in experiment. D) Simulation of spectra for the single layer model 
depicted in E) with uniform decrease of oxide intensities. E) Single layer model with a homogeneous 
(Sii+, Si2+, Si3+, Si4+) surface layer. 



(CA) measurements shows that - within the exper- 
imental error - the roughness and chemical homo- 
geneity of the surface of a native oxide layer on 
Si wafers does not differ from the properties of a 
thick Si02 layer (cf. Table[T]). The analysis of high 
energy X-ray reflectivity measurements on native 
and thick oxide layers leads to electron densities 
that are again similar (Figure |4]). Especially the 
topmost part (0 — 5 A) of the native oxide layer re- 
sembles the density of the thicker layer. 

The stoichiometry of the native oxide layer 
was studied by X-ray photoelectron spectroscopy 
(XPS) combined with Ar ion etching in order to 
reveal the depth profiles for different oxidation 
states/valencies Si'^^ (k — 0, ...,4). Figure [5]a 
shows the Si-2p spectra recorded in normal emis- 
sion mode (take-off angle 0° along the surface nor- 



mal) with two components representing Si and 
Si^+ at lower and higher binding energy, respec- 
tively. Contributions from Si^+, Si^+ and Si^^ 
could not be resolved, but have to be treated as 
a third peak in the background between the Si*^ 
and Si^+ signals. For stepwise ablation of the 
surface, the intensity of the Si^+ peak decreases 
asymmetrically, i.e., it is shifted towards the Si*^ 
signal, forming a shoulder in the intermediate state 
before vanishing (for calibration of ablation see, 
e.g., Ref.^^). For ablation of about 1.5 nm, con- 
tributions from oxide species can no longer be 
observed, which is in accordance with the thick- 
ness of the native oxide measured by other meth- 
ods (Figure |4] and Table [T]). The asymmetry in 
the Si oxide related part of the spectra is charac- 
teristic for intensity contributions from Si^+ and 
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(Si^+, Si^+, Si-^+) when distributed in a double 
layer model as depicted in Figure |5]B. For this type 
of stacking, the simulation of intensity distribu- 
tions within the Si-2p spectra in Figure [5]C shows 
the same characteristics, namely the asymmetric 
decrease of the intensity from the oxide species 
as the experimental data in Figure |5]A. A simi- 
lar asymmetry was reported in previous studies on 
thermally grown oxide layers. For comparison. 
Figure |5]E shows a second model with a homo- 
geneous distribution of the oxidation states within 
the oxide layer. For this scenario, the simulation 
in Figure [5]D predicts a uniform disappearance of 
the oxide contributions, in contrast to the experi- 
mental observations in Figure |5]A. 

In summary, these results show that native sili- 
con oxide layers can be approximated by a dou- 
ble layer system as depicted in Figure [3]A and 
Figure [5]B, where adjacent to the bulk Si a tran- 
sition layer and then a bulk-like Si02 layer fol- 
lows. Since all characterized material properties 
of the Si02 part are similar to the properties of 
thick oxide layers, there are no hints to assume 
different optical properties. The transition layer, 
however, will mos t likely display increased po- 
larizability.'^^IllM] jj^^j means, differences in the 

vdW potentials wafers featuring native and thick 
oxide layers arouse indeed from subsurface contri- 
butions, since the uppermost material is the same, 
namely Si02. Hence, silicon wafers with native 
oxide layers are indeed a good model substrate to 
study van der Waals interactions with coated ma- 
terials. 



Appendix 

Calculation of Hamaker constants 

Based on the Lifshitz theory, the Hamaker con- 
stant can be calculated by"^ 

oo 

A12-32 = -kT^,' ^lMn)^3Mn) (7) 
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Using the Ninham-Parsegian approximation it 
is possible to obtain the from the adsorp- 
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where fj is the strength of an oscillator, cOj its re- 
laxation frequency, and N is the total number of os- 



cillators. For dielectric materials Eq. (10) reduces 

t o 70.71, 

with (Orot the molecular rotational relaxation fre- 
quency (typically in the IR regime), (O^ the main 
electronic absorption frequency (typically ~ 3 ■ 
10^^ s^^ in the UV regime), and £(C0vis) = 'Jyis the 
refractive index in the visible regime. As usually 



(^1 >> tt^-ot (cf. Eq. ^) the first term in Eq. ([12]) 
can be neglected and Eq. (|7]) may be approximated 
by Eq. ([3]) if the adsorption frequencies of the three 
materials are similar. ^ 
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